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By (a) using an expression for the lattice potential of QCD in terms of a continuum running coupling and (b) 
globally parameterizing this coupling to interpolate between 2- (or higher-) loop QCD in the UV and the flux 
tube prediction in the IR, we can perfectly fit lattice data for the potential down to one lattice spacing and at 
the same time extract the running coupling to high precision. This allows us to quantitatively check the accuracy 
of 2-loop evolution, compare with the Lepage-Mackenzie estimate of the coupling extracted from the plaquette, 
and determine the scale ro ten times more accurately than previously possible. For pure SU(3) we find that the 
coupling scales on the percent level for ft > 6. 
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1. Introduction 

The string picture predicts that the static po- 
tential of QCD behaves for large distances like 



V(r) = ar — — 



const 



(1) 



(ignoring string breaking in the presence of dy- 
namical fcrmions), where a is the string tension 
and the "IR charge" e m = 7r/12 for the sim- 
plest string. In the UV V(r) is Coulomb like, 
with a running charge that is known to 2 loops in 
terms of the A-parameter. The only known gen- 
eral method of obtaining quantitative information 
on V(r) is through Monte Carlo (MC) simula- 
tions, which nowadays reach relative accuracies 
of almost 1CU 4 , 1CU 3 at short, respectively, long 
distances. In contrast, the methods used so far 
extract a from these data with an error of several 
percent, A with a much larger error. To obtain 
a one basically fits to an Ansatz of the form ([[]) . 
For A one does something different, namely, one 
defines a running coupling in the force scheme via 
r 2 V'(r) = C F a F (r) (with C F = (N 2 - 1)/2N for 
SU(iV)), and then fits numerical derivatives of the 
potential data to a 2-loop formula for ap(r). 2- 
loop evolution is just becoming good at the short- 
est distances presently available, where however 
lattice artifacts, that one does not know how to 
take into account properly, prevent one from see- 
ing it. With these methods one would therefore 
hardly benefit from more precise data or smaller 



lattice spacings. Here we present a method that 
allows for a unified fit of the whole r range by 
incorporating lattice artifacts and the running of 
the coupling in a fundamental way. Details and 
a complete set of references can be found in Jl| . 

2. The Method 

The general idea is to express the lattice po- 
tential in terms of a continuum running coupling; 
more precisely there are three ingredients: 
• Use the V-scheme, where the running coupling 
is defined via the continuum static potential 



V(q) 



a v (q) 



(2) 



• Parameterize ay (<?) to take into account 2-loop 
QCD in the UV in terms of A, the string pre- 
diction (0) in the IR in terms of a, and to have 
another parameter for the crossover from UV to 
IR. 

• Express the lattice potential as 



V a i(r) = V ~AirC F 



■*/« (2tt) 3 q 2 



(3) 



where aqi = 2sin(a<7i/2) for the usual (unim- 
proved) gluon action. 

The three parameters in ay{q) and the con- 
stant Vo are then fitted by matching (^) to the 
MC data. Though very plausible, (||) must be 
considered an Ansatz. The spectacular success of 
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our fits shows that it incorporates the lattice ef- 
fects at an astonishing accuracy. But before pre- 
senting these fits we describe the 

3. Parameterization of ay (q) 

The following Ansatz for ay(q) satisfies the 2- 
loop /3-function equation and has no Landau pole 
(as long as cq>1,ci>1 and c>0): 



1 



/3 a v (q) 



= In 



l + ^ln fc ( Co + ^A(g)) 



(4) 



where X(q) = ln b (ci + cjp), and b — (3\/ (i^ is a 
ratio of the first two coefficients of the /3-function. 
Of the three dimensionless parameters co,ci and 
c, the first two are fixed in terms of other param- 
eters by matching @) and (Q) to (|l|), 



, b C F K 2 b / 2/3 



1 6+1 

cp hi c 
26 



(5) 



while c is the crossover parameter to be fitted. 

One can show that by iterating the log's in (Q) 
in a suitable way, it is possible, in principle, to in- 
corporate QCD to any number of loops, in terms 
of the (unknown) higher coefficients of the (3- 
function and other, non-pcrturbative parameters. 

4. Results 

The fitting and error analysis is somewhat in- 
volved 0. We here only mention this: Since 
c is quite strongly correlated with A, it is bet- 
ter to use av(q*) instead of A as independent 
fit parameter. Here q* is some UV scale; we 
chose q* = 3.4018/a for easy comparison with 
the Lepage-Mackenzie estimate || of ay{q*) ex- 
tracted from the plaquette Wn. We performed 
fits for various theories, with gauge group SU(3) 
and SU(2), without and with fermions (in the lat- 
ter cases there was no sign of string breaking) . In 
all cases we could obtain x 2 /-/Vdf ~ 1 including 
all points down to r/a = 1, which is impossible to 
achieve with Coulomb + Linear (C+L) type fits. 
In fig. 1 we show one of our fits. For our detailed 
results we again refer to in table 1 we com- 
pare just two quantities with previous estimates 
(we should mention that the /3 = 6.8 data used 
are preliminary H). One is ay(q*), the other the 
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Figure 1. Our fit and MC data for (3 — 6 pure 
SU(3) in the small r region. 



scale ro {§ defined by Cfo;_f(? , o) = 1.65. [The rhs 
of this equation is chosen so that r « 0.5 fm.] 
We can use the exact relation between ay (q) and 
a_p(r) to calculate r . We thereby retain the con- 
ceptual advantages of using a scale like r without 
the errors from derivatives of lattice data. 

In fig. 2 we show our running coupling and com- 
pare it with estimates of ay(q) at various scales 
extracted by the Lepage-Mackenzie method from 
various Wilson loops (the rightmost circle cor- 
responding to Wii) and Creutz ratios. We also 
show the 2-loop approximation to our ay(q). 

5. Discussion and Conclusion 

The good agreement between our more precise 
and the simpler Lepage-Mackenzie estimate of 
ay(q*) increases the confidence in both methods. 

Our fits show at present hardly any sign of sys- 
tematic errors. The only way to check for such 
errors, then, is to come up with other parame- 
terizations of ay{q) that give similarly good or 
better fits. We have done so by taking 3-loop ef- 
fects into account. It seems that slight systematic 
errors exist only at the "edges" of ay(q), i.e. for 
the parameters ay(q*) and a. Those of the for- 
mer are quoted as the second error in table 1. 
In the intermediate region ay(q), and therefore 
ro, does not seem to have significant systematic 
errors. 
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Table 1 



Our and previous determinations of ay(q*) and r 


Group 


n f 


(3 


ay{q*) from 




ro/a from 




our fit 


Wu 


our fit 


C+L fit 


SU(3) 





6.0 


0.1466(19)(8) 


0.1519(35) 


5.30(1) 


5.44(26) 


SU(3) 





6.4 


0.1293(7)(7) 


0.1302(22) 


9.54(3) 


9.90(54) 


SU(3) 





6.8 


0.1140(2)(4) 


0.1153(15) 


14.56(5) 




SU(2) 





2.85 


0.1687(8)(9) 


0.1712(50) 


19.0(3) 


20.6(14)* 


SU(3) 


2W 


5.3 


0.1902(19)(12) 


0.1991(79) 


3.62(1) 


3.7(2) 


SU(3) 


2S 


5.6 


0.1685(14)(9) 


0.1788(57) 


5.20(1) 


5.2(2) 



The av{q*) from Wu is quoted with error ay(q*) i . W, S indicates Wilson, staggered fermions. 
* From simulations using the Schrodinger functional coupling asp(q), see ref. ||. 
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Figure 2. Fit result for ay(q) for (3 — 6 pure 
SU(3) (the solid lines delineating the error), its 
2-loop approximation (dotted), and results from 
the method of ref. @ (circles). a- 1 = 2.0 GeV. 



In contrast to yfo j A, the quantity roy/a scales 
very well; it equals about 1.14 for the nf — 2 theo- 
ries of table 1, 1.17 for pure SU(2), and 1.17-1.19 
for pure SU(3). For the latter tq^/g therefore 
scales at the 1% level for j3 > 6.0. At the edges 
of a,y(q) the scaling violations are slightly larger, 
but we expect || more precise data and analysis 
to further decrease them. 

As one might suspect from fig. 2 and table 1, 
we are beginning to gain control of ay (q) at the 
fraction of percent level on all momentum scales 
(except perhaps in the far IR). For intermediate 



momenta below about 5 GeV this is basically al- 
ready the case; for larger momenta 2-loop evolu- 
tion is becoming good, so there is hope to soon 
achieve this accuracy also in the UV. 

The bad scaling properties of quantities involv- 
ing A are not surprising: We found that the error 
of A is quite large — without leading to a large 
error in ay (q) itself — because A is strongly cor- 
related with higher order and non-perturbative 
contributions at an intermediate range, appear- 
ing here in the form of the fit parameter c. 

For the future we note that our method can also 
be used to fit results from improved gluon actions 
(with the obvious modifications in eq. (§)), and 
can be extended to incorporate string breaking. 

Finally, it should be clear that our ay(q) is 
useful for potential models of heavy quarks and 
anywhere one wants to extend perturbation the- 
ory without running into the Landau pole. 
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